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Abstract 

Duality relations for the 2D nonhomogeneous Ising model on the finite square lattice 
wrapped on the torus are obtained. The partition function of the model on the dual lattice 
with arbitrary combinations of the periodical and antiperiodical boundary conditions 
along the cycles of the torus is expressed through some specific combination of the partition 
functions of the model on the original lattice with corresponding boundary conditions. It 
is shown that the structure of the duality relations is connected with the topological 
peculiarities of the dual transformation of the model on the torus. 
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The duality relation for the two-dimensional Ising model was discovered by Kramers and 
Wannier [1]. In their work, Kramers and Wannier showed the correspondence between the 
partition function of the model in low-temperature phase and the partition function of the 
model on the dual lattice in high-temperature phase and vice versa: 

(cosh2K))- N Z{K) = (cosh2K)- N Z{K) (1) 

sinh2K • sinh2K = 1. 

Using this self-duality property, the critical temperature of the 2D Ising model was determined 
[1] before Onsager had obtained the exact solution [2]. 

Kadanoff and Ceva [3] generalized the Kramers- Wannier duality relation (1) to the nonho- 
mogeneous case (the coupling constants are arbitrary functions of lattice site coordinates) with 
spherical boundary conditions 

Y[(cosh2K,(r)y 1/2 Z[K] = ]J(cosh 2 K, (r)y l /2 Z[K], (2) 

smh2K x (r) • smh2K y (r) = 1, smh2K y (r) ■ sinh2K x (r) = 1. (3) 

Here /i — x,y and r, f, K^r), K^if) are coordinates and coupling constants on the original and 
dual lattices respectively. Since the Kadanoff-Ceva relation (2) defines the connection between 
functionals, this relation is very useful for analysis of the thermodymamic phases of the model. 
Thus, for example, this relation allows one to define correctly the disorder parameter //, to 
obtain the duality relation connecting correlation functions on the original and dual lattices, 
to define "mixed" correlation functions (crfa) . . . a(rj)fj,(rk) ■ ■ ■ M r 0) an d so on (see Ref. [3]). 

As was already mentioned in Ref. [1,3], relations (1) and (2) can not be understood literally. 
So, for example, using the method of comparing high- and low-temperature expansions for 
deriving the duality relation (1) in the case of the periodical boundary conditions, it is hard to 
take into account and to compare the graphs wrapping up the torus. In fact Eq. (1) is correct 
in the thermodynamic limit (for the specific free energy). However for the nonhomogeneous 
case the procedure of thermodynamic limit is rather ambiguous. In Ref. [3] the duality relation 
(3) was obtained for spherical (nonphysical for the square lattice) boundary conditions. 

Since duality is a popular method of nonperturbative investigation in field theory and sta- 
tistical mechanics (for review see Ref. [4]), it is important to formulate a duality transformation 
for finite systems. Recently, we have suggested [5,6] the exact duality relations for the nonho- 
mogeneous Ising model on a finite square lattice of size N = n x m wrapped on the torus: 

n(cosh2i^(f)r 1/2 Z[K] = n(cosh2i^(r))- 1/2 T Z[K\. (4) 
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and components of the four- vector (Z[K] for the dual lattice) 

Z[K] = (^ZM.ZM^M), 

are partition functions Z^ a '^[K] (a,/3 = a,p) of the Ising model with corresponding combina- 
tions of the periodical (p) and antiperiodical (a) boundary conditions along the horizontal X 
and vertical Y axes: 

zW[K\=Y,<r PH ™ [KM , (6) 



- /Hf^tf, a] = £ a(r)(K x (r)V a x + 7^(r)V>(r), 



(7) 



where r = (x, y) denotes the site coordinates on the square lattice of size N — nxm, x — 1, . . . ,n 
y = 1, . . . , m; c(r) = ±1; K x (r) and K y (r) are the coupling constants along the horisontal X 
and vertical F axes respectively. The one-step shift operators V x , V y act on a(r) in the 
following way 

V x a(r) = a(r + x), V y a(r) = a(r + y), 

where x, y are the unit vectors along the horisontal and vertical axes. For the periodical 
(antiperiodical) boundary conditions along X and Y axes we have 

V^ (a) a(n, y) = +(-)<r(l, y), Vf m) = +(-)a(x, 1). 

We denote site coordinates, functions and functionals on the dual lattice by "tilda" : f, cf(f), 
K^(r), H[K,a], Z[K], .... A site coordinate on the dual lattice coincides with a coordinate 
of the plaquet center on the original lattice: 

r = r + (x + y)/2. 



(8) 



In Ref. [5] the duality relation (4) was proved for homogeneous and weakly nonhomogeneous 
distributions of the coupling constants. We also have checked the duality relation (4) for lattices 
of small sizes by direct calculation on the computer. As a corollary of Eq. (4), we obtained 
[5,6] the duality relations for the two-point correlation function on the torus, for the partition 
functions of the 2D Ising model with magnetic fields applied to the boundaries and the 2D Ising 
model with free, fixed and mixed boundary conditions. 
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In this Letter we would like to propose a simple way of the derivation of the duality relation 
(4). For the derivation it is convenient to use the representation of Hamiltonian (7) in terms 
of Hamiltonian H^'^ of the Ising model with the magnetic dislocation corresponding to the 
boundary conditions (a, (3) and the periodical boundary conditions for V M : 

r,n 

Here the coupling constants configurations [lf( a '^] define the following magnetic dislocations: 
(i) For the Hamiltonian H^'^: 



KM(r>) = K x (r>), 

K M(r') = -K y (r'), if r'ef; 

KM(r) = K,(r), if r^B^; 



where we introduced the denotion for the following boundary cycles on the torus 

B x = {{x,i),x = l,...,n}, i = l,m, 

(ii) For the Hamiltonian Hp' p ^\ 

( K^(r') = -K x (r'), 
D Y = I Kjw)(r>) = Ky (r'), if r' e Bp>; 
[ Kjw)( r ) = K^r), ifr££^ n) ; 

where the path By denotes the other boundary cycles on the torus 

By = {{i,y),y = l,--,m} , i = l,n, 

(iii) For the Hamiltonian Hp' a ^ we have Dxy = Dx U Dy. 

It is evident that H^' p ^ has not magnetic dislocation. Nevertheless, let us denote configu- 
ration of the coupling constants in this case as -Do and call coupling constants configurations 
Dq, Dx, Dy, Dx,y as the basic magnetic dislocations for corresponding H^'^K 

Note that the partition function (6) is invariant with respect to the Z 2 local gauge trans- 
formation U[t] [7,8]: 

U[t] [K,a]= [K',a']= \{K,(r) = r(r)K^r)r(r + £)}, {a\r) = r(r)a(r)}}, (9) 
where r(r) = ±1. Let us apply arbitrary gauge transformation to the partition function Z^ a '^ 
Ulr^eM-pH^lK, a}} = ^expf-^^fX', a'}} = Z (a ^[K']= Z (a ^[K]. 
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Here new coupling constants configuration [K'\ can contain both the deformation of the cor- 
responding basic magnetic dislocations and new closed magnetic dislocations. Let us denote 
by Qb' a \ [}( a ' p \ il^'^ the classes of gauge equivalent configurations [K'\ generated by the 

gauge transformations from the basic magnetic dislocations D , Dx, Dy, Dx,y respectively. It is 
evident that these classes have not intersections because the basic magnetic dislocations are the 
homotopy-nonequivalent paths on the two-dimensional torus (two arbitrary coupling constants 
configurations from different classes QS a ^ can't be connected by the Z 2 gauge transformations). 

For the dual transformation of the partition function (6) let us use the standard method [4] 
of the passage to dual spins. At the beginning we consider the dual transformation of Z^ p ^\ 

z (p,p) [k] = X)expE^(r)^(r)V^(r)] = 

[a] r,n 

EniE^w(^M)K^H/i)) Wr) ]= 

E E II fl,(r) (K, (r)) U(^)) m = ^ E II P Ur)(K,(r)) J[ (10) 
M H r >» r [i„] r >v r 

where 

P Wr) (i^(r)) = cosh(^(r))exp(/ M (r)lntanh^(r)) (11) 

and 

^(r) = l x (r) + l y (r) + l x (r — x) + l y {r - y), l^r) = 0, 1. (12) 

Here 82(1^) is a Kronecker 5-function mod 2 : it is zero if ip is odd and one if ip is even. Note 
that in Eq. (10) we have a product of the (^-functions with linking arguments (the same 1^ 
is contained in two 5-functions) . In order to solve constraints generated by the product of 
(^-functions £ M (r) is usually expressed through dual spins a(r) [4] 

i M (r) = ^(l-ar(r)ar(r-i/)), p±v. (13) 

Substituting Eq. (12) in Eq. (9) we obtain 

Z M[K\ = l^n( sinh2 ^W) | expE5(f)^(f)V^(f)] = 

^(cosh2K u (r))- 1 Z^[K], (14) 

where the dual coupling constants K v (r) are defined by (3) or 

tanhi^(r) = e~ 2 ^\ n^v. (15) 
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To derive Eq. (14) we used identity 

cosh 2 2K fl (r) cosh 2 2i^(f) 
sinh2i^(r) ~ S mh2K u (r) ' 



Since in Eq. (14) we sum over [a] (this counts each configuration [l^] twice because [a] — > [—a] 
gives the same [1^]), we must introduce factor 1/2. 

However it is not hard to note that we can construct many other solutions for for which 
the coupling constant configurations are connected with [K] in (14) by means of the gauge 
transformation U[f] on the dual lattice (see Eq. (9)). Such configurations form a class Q p,p 
of gauge-equivalent coupling constant configurations . By analogy with the Ising model on 
the original lattice wrapped on the torus we must expect the existence of solutions for Z M 
which lead to homotopy-nonequivalent classes VL^ P,P \ Ct^ p,a \ Q( a ' p \ Cl( a ' a ' ) of dual coupling con- 
stants configurations. Really it is easy to write the solutions for which lead to the basic 
magnetic dislocations D x , D Y , D X y in Vt^ p,a \Vt <ya,p \Vt <ya,a "> respectively. So, taking into account 
Eq. (8), we have for D X - 

l x (r) = ±(l + a(r)a(r-y)), (16) 
l y (r) = l-(l-a(r)a(r-x)), if rG 5?; 
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l li (r) = ±(l-a(r)a(r-V)), iir^B x ] -, 



for D Y : 



and for D X y\ 



l v {r) = ^{l + a(r)a(r-x)) 1 (17) 
l x (r) = ^(l-a(r)a(r-y)), if rG B?; 
Ur) = l(l-a(r)a(r-v)), iir^B?; 



l x (r) = ±(l + <r(r)a(r-y)), (18) 
l y (r) = ±(l-a(r)a(r-x)), if rG 5§>; 

l v (r) = ±(l + a(r)a(r-x)), 
l x (r) = ±(l-a(r)a(r-y)), iireB?; 
W) = 1(1 - a(r)a(r-u)), if r # B x \b£\ 



In Eqs. (16)-(18) n^v. 
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Since arbitrary solutions satisfying the product of ^-functions in Eq. (10) lead to the dual 
coupling constant configurations which contain the finite number of closed magnetic dislocations 
on the dual lattice, it is obvious that these configurations can be generated by the gauge 
transformations from the basic magnetic dislocations and occur in the corresponding class 
Q(a,P) m Then in Eqs. (10), (14) at transformation to the dual lattice we must sum over the 
homotopy-nonequivalent solutions (13), (16)-(18). Taking into account that basic magnetic 
dislocations D , Dx, Dy, Dx,y lead to the partition functions with corresponding boundary 
conditions, we obtain: 

Y[(cosh2K tl (r))-^Z^[K] = 

r,fi 

- H(cosh2K u (r))~^ (z M [K] + Z^ a \K\ + Z {a ^[K] + Z^[K]) . (19) 
_ — 

Now, let us consider the dual transformation of partition function Z^ p,a \ As we have discussed 
above, it contains magnetic dislocation D X - Again, doing the duality transformation, it is 
necessary to sum over the homotopy-nonequivalent solutions (13), (16)-(18). As a result we 
obtain the duality relations similar to Eq. (19) but with the minus signs before Z^ a \K\ and 
Z( a > a ) [K] . Appearence of these signs is connected with the presence of magnetic dislocation 
D x in Z^ v > a \ Let us show this. Considering the transformation to the dual lattice in the same 
way as in Eqs. (10), (14), we obtain the partition function Z^ p,p \K^ p,,p) ] with the following 
dislocation (see Ref. [3]): 

f KjP*)(?))=K x (r , )+q, 
G x = I K^(r') = K y (r'), if r> G 

[ K^(r)=K ll (r), iff ? B™; 

where the path B^ denotes the following boundary cycle on the dual torus 

B { x ] = {(x,m),x = 1, ...,n} , 

that is the dual transformation transforms dislocation Dx to dislocation Gx- 

The partition function Z^ [K^] can be written through the following correlation function 

[3] 

i n < a(l,m)a(2,m)a(2,m)a(3,m)...a(n,m)a(n + l,m) > Z {p ' p) [K} = i n Z M [K], 

where Z^ P ' P \K\ does not contain magnetic dislocations. In other hand, observing that except the 
dislocation G x the partition functions Z (j ' ,a ^[K^ p,a ^\ and Z ( - a,a ^[K ( - a ' a ' > } have magnetic dislocation 
D y , we obtain, for example, for the last partition function 

i n < d(l,m)a(2,m)a(2,m)a(3,m)...a(n,m)a(n+ l,m) > Z M [K] = -i n Z {a ' a) [K}. 
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Here the minus sign appeares because of the antiperiodical boundary conditions. Similarly the 
minus sign appears before Z^' a ^[K]. Then, taking into account these signs, one gets 

n(cosh2i^ ^r)y^Z M [K] = 

r,n 

- ;Q(cosh2^(f)H (Z {P ' P) [K] - Z M [K] + Z ia ^[K] - Z^[K]) , (20) 

2 — 

r,u 

where the factor i n is cancelled by the same factor appearing from the normalizing product 
T7 cosh 1 / 2 K. The same way, considering the dual transformation for Z 1 -" 1 ^ and Z^'^, we obtain 

n(cosh27^(r)r 1/2 Z[K] = n(cosh2i^(f)r 1/2 f Z[K}. (21) 

Multiplying the right and the left sides of Eq. (21) by matrix T (5), one obtains Eq. (4). 

Thus the structure of duality relation (4) is connected with the topological peculiarities 
of the dual transformation of the Ising model on the torus. From here, one can assume that 
dual relations for other models with Zi symmetry, for example, the eight-vertex model, have 
the similar structure. The method of derivation of duality relation (4) suggested in this paper 
can be generalized for the two dimensional lattice models with Z N symmetry. Results of this 
reseach will be published in the following paper. 

The authors are indebted to A. Belavin and A. Morozov for continuous support and are 
grateful to M. Lashkevich and A. Mironov for useful discussions. One of the authors (V.S.) 
thanks A. Morozov for the hospitality and the exellent conditions at ITEP where this work has 
been begun. This work is partially supported by grant INTAS-93-1038. 
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